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Abstract 

We explore the discovery potential of cosmic ray physics experi- 
ments for Standard Model processes involving the nonperturba- 
tive production of > 0{a^) ~ 30 weak gauge bosons. We 
demonstrate an experimental insensitivity to proton-induced pro- 
cesses and emphasize the importance of neutrino-induced pro- 
cesses. The Fly's Eye currently constrains the largest region of 
parameter space characterizing multi-W phenomena if a cosmic 
neutrino flux exists at levels suggested by recent models of active 
galactic nuclei. MACRO (DUMAND) can constrain or observe 
additional regions by searching for 1-100 (1-10) characteristic 
near-vertical (near-horizontal) spatially compact energetic muon 
bundles per year. 
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1 Introduction 



Lowest-order perturbative calculations in the Standard Model for pro- 
cesses involving the production of > 0{a^) ~ 30 weak gauge bosons suggest 
an explosive (and unitarity violating) growth of the associated parton-parton 
cross section above center-of-mass energies > 0{a^Mw) — 2.4 TeV [1-6]. 
Not surprisingly, this intriguing circumstance has drawn considerable interest 
in attempts to either substantiate or dismiss the possibility that the lowest- 
order result is a harbinger of spectacularly rich phenomena which may be 
observable at the next generation of hadron colliders (see ref. [7] for an 
overview). When unitarity is restored (most likely due to nonperturbative 
effects) it is conceivable that the cross section for multi-W phenomena^ is un- 
observably small. However, despite concentrated theoretical efforts toward 
resolving the issue of large versus small cross sections no definitive answer 
has arisen; there is a very real possibility that the question may first be set- 
tled experimentally. Given the stakes involved, a quantitative consideration 
of experimental constraints on multi-W production is clearly desirable. 

The high-energy, high-luminosity environment of hadron colliders such as 
the proposed Large Hadron Collider (LHC) and the Superconducting Super 
Collider (SSC) would be ideal for observing or constraining multi-W phe- 
nomena over a wide range of energies and cross sections [T^. However, 
until such machines are commissioned, cosmic rays provide the only access to 
the required energy scales. In this paper we adopt a purely phenomenological 
approach and explore the feasibility of discovering or constraining multi-W 
phenomena in the context of cosmic rays. We consider both atmospheric 
and underground phenomena induced by cosmic ray protons and neutrinos. 
Some of the consequences of neutrino-induced phenomena for future under- 



water detectors such as DUMAND W^ and NESTOR have already been 



^It is understood in the following that by 'multi-W phenomena we refer to processes 
producing a total of > 0{a^) ~ 30 W's and Z's. 
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discussed in refs. |T3|, [T^. 

By adopting a simple parameterization of multi-W phenomena on the 
parton level, we evaluate the discovery potential of various experimental ar- 
rangements over a space of possible theories involving multi-W production. 
We demonstrate that even if proton-induced atmospheric multi-W phenom- 
ena occur in Nature, the features of the resulting air showers are unlikely 
to allow one to distinguish them from fluctuations in a much larger back- 
ground of generic showers. On the other hand, ultrahigh energy neutrinos, 
for which a sizeable flux has recently been conjectured from sources such 
as active galactic nuclei [15-18], offer exciting possibilities for observing or 
constraining multi-W phenomena. Subsurface detectors such as AMANDA 
1|], DUMAND, MACRO |g, NESTOR, and NT-200 m can be sensitive to 



neutrino-induced multi-W phenomena in some regions of multi-W parameter 
space. 

The structure of the paper is as follows. In sect. 2 we characterize multi- 
W phenomena by a two-parameter working hypothesis which frees us from 
specifying an underlying (most likely nonperturbative) mechanism for multi- 
W production. We also describe the gross features of multi-W phenomena 
and present discovery limits for the LHC and SSC. We discuss proton-induced 
and neutrino-induced multi-W phenomena in sects. 3 and 4, respectively. 
In each section we consider a variety of detection techniques and present 
discovery limits which may be contrasted with the superior sensitivity of 
future hadron colliders. In sect. 5 we summarize our results and conclude. 
In an appendix we outline our quantitative description of multi-W processes. 



2 Working Hypothesis 

In the absence of a reliable first-principles calculation of multi-W pro- 
duction, we adopt a working hypothesis which allows us to parametrize the 
essential features of nonperturbative phenomena without committing our- 
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selves to a specific underlying mechanism. We model the parton-parton cross 
section for multi-W production by 

aoe{V~s-Jio). (1) 



0" multi- W 

For the purposes of this paper a parton is any weakly interacting particle 
(for example, q,i^,e). The 9— function in eq. |l| embodies an idealized onset of 
nonperturbative multi-W production above a parton-parton center-of-mass 
threshold energy of \/so with a cross section of (Tq. While these definitions are 
convenient, it should be kept in mind that the true behaviour of the parton- 
parton multi-W cross section near threshold may be more complicated than 
eq. IH; in that case y/so and ctq in eq. ^ should be interpreted as effective 
parameters. 

We wish to explore the possibility of observing multi-W phenomena char- 
acterized by parameters in the range 



~ 2.4 TeV < Vso < 40 TeV, 
aw 

(2) 

4^^100pb < do < crZixC-^y ^lOfih. 

The lower limit of a/Iq is suggested by the energy scale at which perturbation 
theory becomes unreliable [§, ^, H whereas the upper range is of the order 
of the sphaleron mass p2[ . The lower range of (Jq follows from dimensional 



arguments, being characteristic of a geometrical "weak" cross section. The 
upper range of (Tq is a geometrical "strong" cross section suggested by analo- 
gies between the weak SU(2) gauge sector and the color SU(3) gauge sector 
[|]. Admittedly, our current theoretical understanding of weak interactions 
renders ctq > 1 fih an unlikely scenario; we nevertheless include it in our 
analysis so that it may be tested experimentally. For definiteness, we will 
assume throughout this paper that ctq refers to the production of exactly 
30 W bosons; allowing for the production of variable numbers of W's (and 
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Z's and possibly prompt photons) is straightforward but is an unnecessary 
comphcation at the level of our investigation. 

Figure ^ shows the multi-W production cross section for protons and 
neutrinos (energy E) striking stationary nucleons and electrons. The curves 
are obtained by convoluting eq. |l] with the corresponding quark distribution 
functions inside nucleons (see appendix for details). The results are universal 
in the sense that they have been scaled by o"o and ^/Iq- 

The simultaneous production of (9(30) W bosons at future hadron collid- 
ers like LHC or SSC would lead to spectacular signatures [9,10]. Since the 
average number of charged hadrons (mainly vr^'s) from hadronic W decays 
is {n^J^^^^'^^"^^^) ~ 20, one could typically expect 

30 X Br{W ^ hadrons) x (ni^^''"''™"'^) ~ 400 vr^ (3) 

in one multi-W event accompanied by ~ 400 photons from the decay of 
^ 200 vr^'s. The charged hadrons would have a minimum average transverse 
momentum of order 

> O{mw/S0) ~ (2 - 3) GeV, (4) 

if the W bosons are produced without transverse momentum. Similarly, one 
could expect ~ 5 prompt muons (~ 3 from W decays and ~ 2 from c, b, or 
r decay) carrying a minimum average transverse momentum of 

> 0{mw/2) - 40 GeV. (5) 

A similar situation holds for other prompt leptons such as e^,^ etc. It is 
hard to imagine that any other process in the Standard Model can mimic 



such a final state IG 



Figure ^ shows the regions in a/Iq — ctq space accessible to the LHC 
(Vs = 14.6 TeV; C = 10^^ cm'^s'^) and the SSC (^s = 40 TeV; C = 
10^^ cm~^ s~^). The contours correspond to 1 and 10 events (assuming 100% 
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detection efficiency) for 10^ s of operation. These contours may be used 
as a benchmark to evaluate the effectiveness of various cosmic ray physics 
experiments for constraining muhi-W phenomena. 

3 Proton-Induced Multi-W Processes 



Cosmic ray protons and heavy nuclei constitute a guaranteed flux of 
high-energy primaries potentially capable of initiating multi-W phenomena. 
In this section we explore the possibility of exploiting this cosmic flux and 
isolating multi-W phenomena from generic hadronic reactions. We restrict 
our attention to cosmic ray protons since they provide the dominant flux in 
terms of energy per nucleon. 

Since the flux of proton-induced multi-W air showers is anticipated to be 
small for even the most optimistic scenarios, one must resort to experiments 
with large effective areas and/or long exposure times. We will discuss three 
types of relevant experiments: 

1) conventional surface arrays which measure the e,7,/i, etc., content of 
air showers (AGASA [^, CYGNUS GAS A |2§, EAS-TOP EAS- 
100 13, HEGRA [ID, KASGADE H], etc.), 

2) underground experiments sensitive to downward through-going TeV 
muons (LVD fSC 



MAGRO 20 



Soudan-2 etc.) or underwater (-ice) 
experiments which detect energetic muons by Gherenkov light (AMANDA 
DUMAND ig, NESTOR NT-200 |2l|], etc.), and 



3) the Fly's Eye an optical array which is sensitive to nitrogen fluo- 
rescence light from air showers. 



3.1 Conventional surface arrays 

Using a network of detector elements, conventional surface arrays recon- 
struct the features of an extensive air shower by interpolating or extrapolating 
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measurements of a shower's particle content. Among the particles sampled 
are photons and electrons with Eg, > 0{1 MeV), muons with Efj_ > 0{1 
GeV), and, in some experiments, hadrons with -Ehad ^ 1 GeV. 

Consider a surface array (area A) sensitive to showers above a threshold 
energy -Ethresh- For proton- induced phenomena it follows from elementary 
considerations that the number of multi-W showers occurring during time t 
is 

Number of multi-W showers = tAQ T dE ^"^"|^tj-w(-^) ^-^^^^ ^ 

^thresh 0"mel(-^) 

where jp{E) = dNp/ (dAdtdQdE) is the cosmic proton flux and Q is the solid 
angle acceptance of the array. For the cases we consider the inelastic proton- 
nucleon cross section af^i is dominated by generic QCD interactions and 
may be approximated by afj^i{E) ~ const. ~ 100 mb. In principle, multi-W 
production through proton-electron collisions is also possible. However as 
illustrated in Fig. |I|, the cosmic proton threshold energy for multi-W produc- 
tion in pe~ collisions is nip/ me — 1800 times larger than the corresponding 
pN threshold. Moreover, the pe~ multi-W cross section is at least 100 times 
smaller than the pN multi-W cross section; hence we neglect multi-W pro- 
duction through pe~ collisions. 

Figure § shows contours for the number of proton-induced multi-W air 
showers at zenith angles 9 < 60° striking a 100 km^ conventional surface 
array in 10^ s. For our calculations we use the cosmic proton flux of the Con- 
stant Mass Composition (CMC) model (see Fig. H). Though the total 
flux of ultrahigh energy hadronic cosmic rays is relatively well measured, its 
composition {i.e., percentage of p, Fe, Mg etc.) is somewhat less certain. Our 
use of the CMC proton flux introduces an inherent, though not critical un- 
certainty in this respect; the CMC flux is a compromise between proton-rich 
and proton-poor |35[ scenarios. For purposes of illustration we optimisti- 
cally assume ii^thresh = 1 PeV which accommodates all multi-W thresholds 
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above > 2.4 TeV. In 100 km^ arrays like AGASA and EAS-100, inter- 
detector spacing on the order of .5-1 km makes -Ethresh = 100 — 1000 PeV 
more realistic but does not change our conclusions. 

Though the region of the {y/so, ao) plane accessible to surface arrays is 
clearly limited by the absolute rate of proton-induced multi-W phenomena, it 
is instructive to consider how one might distinguish an air shower of multi-W 
origin from a generic hadronic air shower. For the remainder of this section, 
we restrict our attention to the optimistic scenario of parton-parton multi-W 
threshold of ^/so = 2.4 TeV with o"o = 10 /ib; for this choice of param- 
eters a 100 km^ array would see approximately 110 (45) multi-W showers 
in 10^ s at zenith angles < 60° for ^thresh = 1 PeV (100 PeV). As il- 
lustrated in Fig. H, the combination of a rising cr^^j^i_w(-^) with a falling 
cosmic proton flux spectrum implies that typical multi-W showers would 
have energies well above the corresponding cosmic proton threshold energy 
of E^'^^'^^ = So/(2mp) ~ 3 PeV. The most probable shower energy is ~ 30 
PeV and the average shower energy is ~ 250 PeV due to a long tail on the 
distribution. 

Consider the characteristics of the most probable ( ~ 30 PeV) multi-W 
air showers. To phrase our results in experimentally relevant terms we use 
the computer program SHOWERSIM to simulate multi-W air showers and 
generate samples of generic proton-induced and iron-induced showers (see 
appendix for details). Figure H compares 30 PeV multi-W, proton and iron 
showers in terms of radial particle densities (with respect to a vertical shower 
axis) of electrons {Eg > 1 MeV), muons (i?^ > 1 GeV) and hadrons (-Ehad > 
1 GeV). Each curve is averaged over 25-100 showers taking into account 
the distribution of the depth of first interaction in the upper atmosphere. 
The densities in Fig. ^ correspond to an an observation depth of 800 g/ cm^ 
(roughly the CYGNUS array depth l2^). The corresponding average particle 
numbers are listed in Table 1. We neglect possible systematic uncertainties in 
the experimental determination of shower energies which may be important 
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in practice. 

The differences between the particle density profiles of 30 PeV showers in 
Fig. ^ are hardly striking. While there are identifiable systematic differences 
between average showers of different origin, the differences do not appear to 
be sufficient to discriminate between multi-W showers and fluctuations in 
generic proton or iron showers. We emphasize this point by noting that in 
the CMC flux model, the differential fluxes of 30 PeV generic proton-induced, 
iron-induced and multi-W showers (with ^/so = 2.4 TeV, o"o = 10 /ib) stand 
in the proportion p : Fe : multi-W ~ 1.2 x 10^ : 1.1 x 10^ : 1. 

Considering the spectacular underlying nature of multi-W phenomenon 
the similarities between proton-induced multi-W showers and generic proton 
and iron may appear surprising. However, a proton-induced multi-W shower 
is actually a superposition of an "interesting" prompt shower component 
seeded by the instantaneous decays of W bosons and a generic "uninterest- 
ing" component initiated by the proton fragment which does not participate 
in the multi-W production subprocess. The proton fragment emerges from 
the region of multi-W production, hadronizes, and subsequently generates a 
generic hadronic shower deeper in the atmosphere. Since the proton fragment 
typically carries a substantial fraction of the primary proton energy, the so- 
called leading particle effect, a large generic component to a proton- induced 
multi-W shower jeopardizes the chances of isolating multi-W showers from 
"pure" generic showers. As shown in Fig. |^, approximately 60% (20%) of 
the total energy in a 30 PeV (5 PeV) proton-induced multi-W air shower is 
carried by the generic component. 

To minimize the effects of the uninteresting generic component of multi- 
W showers and accentuate the prompt component one may be tempted to 
consider multi-W showers close to the relevant threshold energy. Flatter par- 
ticle densities near the core of 5 PeV multi-W showers (compared to 30 PeV 
showers) in Fig. ^ suggest the larger transverse momentum characteristic of 
particles from W decay. However, it is unlikely that one can capitalize on 
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such effects in practice. Due to the slow turn-on of cr^uiti-w (^^^ ^ig. |1|), 
showers just above threshold are rare; for the case at hand only .0002% (17%) 
of proton-induced multi-W showers have energy less than 5 PeV (30 PeV). 
This corresponds to to only one proton-induced multi-W air shower (with 
-^'shower < 5 PcV) in a 100 km^ array every 4400 years! In summary, the 
prospects for detecting proton-induced multi-W phenomena using conven- 
tional surface arrays are poor. 



3.2 Underground detectors 

In an effort to overcome the complications introduced by a large generic 
shower component in proton-induced multi-W air showers, it is helpful to 
concentrate on aspects of air showers which reflect the nature of the primary 
hard interaction. Energetic muons produced from the prompt decays of W 
bosons or from the weak decays of energetic mesons are good candidates 
in this respect. For definiteness, we restrict our attention to muons with 

> 1.5 TeV; such muons can penetrate to deep underground detectors 
such as MACRO |20|. 

As shown in Fig. ^ the lateral distribution of TeV muons in 30 PeV 
multi-W showers is flatter than the corresponding distributions for 30 PeV 
generic p-induced or Fe-induced showers. The flatter multi-W distribution is 
characteristic of the large transverse momentum of the prompt muons from 
W (and Z) decays, eq. ^, and the decays of pions and kaons, eq. |[ The 
corresponding average number of TeV muons in each type of shower is given 
in the last column of Table 1. 

Though distinct, TeV muon signatures from proton- induced multi-W phe- 
nomena are limited by small event rates and practical detector sizes. The 
event rate contours for MACRO can be estimated from the 100 km^ sur- 
face array contours in Fig. ^ simply by scaling by the area ratio (12 m x 
77 m)/100 km^ ~ 10~^, implying negligibly small rates. Even for poorly mo- 



9 



tivated parameters such as (v^ = 2.4 TeV, ctq = 100 ^h) one would need 
an underground detector with an area sensitive to downward moving muons 
of O{10^ m^) to see penetrating muons from one proton- induced multi-W 
event in 10^ s. In view of these small rates, present and future underground 
detectors are not sensitive to penetrating muon signatures of proton-induced 
multi-W phenomena. 



3.3 Fly's Eye 

Finally, we turn to the discovery potential of the Fly's Eye p2 



an 



optical array sensitive to nitrogen fluorescence light from air showers whose 
trajectories do not necessarily intersect the array. By detecting fluorescence 
light emitted as air showers streak across the sky, the Fly's Eye is capable 
of reconstructing the longitudinal development of air showers with energy 
greater than -Ethresh = 100 PeV (see ref. for a pedagogical introduction). 

Analogous to our previous calculation for conventional arrays, the number 
of proton-induced multi-W showers seen in time t by the Fly's Eye is given 
by 

/■oo (tP^ ( E) 

Number of multi-W showers = t dE "^""l'-^^ ' jp{E) An{E) , (7) 

(^) 

where the acceptance, AQ, is a function of energy and hence appears under 
the integral. Fig. ^ shows event number contours for proton-induced multi- 
W air showers corresponding to 10^ s operation of the Fly's Eye using the 
CMC proton flux and the acceptance of ref. Due to its sensitivity, the 
Fly's Eye operates only on clear, moonless nights; approximately 2 x 10^ s 
of observation time is possible in one calendar year |^ . 



Despite the limited region of multi-W parameter space accessible to the 
Fly's Eye through proton-induced multi-W air showers we consider whether 
multi-W showers, if present, can be differentiated from generic showers. In 
the same spirit that one expects generic Fe-induced air showers to develop 



10 



more rapidly than generic p-induced showers, one might expect that the large 
initial multiplicity in a multi-W event (from the immediate decay of the W 
bosons) leads to an accelerated development of the corresponding air shower. 
To test this hypothesis we compare the longitudinal profiles {i.e., the number 
of electrons as a function of shower depth) of multi-W showers with those 
of generic p- and Fe-induced air showers. Fig. ^ shows samples of 150 PeV 
vertical showers of each type. The multi-W air showers assume a parton- 
parton threshold of a/Iq = 5 TeV so that 150 PeV showers are approximately 
a factor of 10 above the corresponding proton threshold of So/(2mp) ~ 13 
PeV and are not atypical; for this choice of threshold the Fly's Eye would 
expect to see 1-10 multi-W air showers within 10'' s observation time for o"o 
in the range 10-100 /ib. 

For ease of comparison in Fig. ^, the depth of first interaction of 150 PeV 
cosmic ray protons (Fe nuclei) was fixed at 42 g/cm^ (11 g/cm^) which corre- 
sponds to the average depth of first interaction. While systematic differences 
are evident between profiles of different origin, the longitudinal profiles of 
multi-W showers are not sufficiently distinctive to prevent confusion with 
fluctuations in generic air showers. For (^/sq = 5 TeV, ctq = 10 /ib) the dif- 
ferential fluxes of 150 PeV generic proton-induced, iron-induced and multi-W 
showers stand in proportion to ~ 80000 : 75000 : 1 in the the CMC model. 

4 Neutrino-Induced Multi-W Processes 

A cosmic flux of ultrahigh energy neutrinos would a provide a novel op- 
portunity to search for multi-W phenomena. Unlike proton-induced multi-W 
production which must compete with (9(100 mb) generic hadronic processes, 
neutrino- induced multi-W production competes only with 0{nh) weak inter- 
action processes. Fortuitously, considerable enthusiasm has been generated 
recently by predictions of a large flux of cosmic neutrinos from active galac- 
tic nuclei (AGN) [15-18]. Indeed, in the model of Stecker et al. (see Fig. ^) 
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the diffuse ffux of PeV neutrinos from AGN is comparable to tlie flux of 
cosmic protons in the CMC model! If only generic charged current interac- 
tions are operative, the Stecker et al. flux suggests that DUMAND should 
see 154(66) single muon events per year with energies E'^ > 100 GeV (10 
TeV) at zenith angles 9 > 70° [l^. In this section we explore alternatives 
for uncovering neutrino-induced multi-W phenomena assuming that a flux of 
ultrahigh energy neutrinos does, in fact, exist. 

We restrict our attention to multi-W phenomena in neutrino-nucleon col- 
lisions and neglect multi-W phenomena from neutrino-electron collisions. As 
may be deduced from Fig. |I|, the ratio c'"muiti- w/'^muiti- w reaches a maximum 
of ~ 1/30 when cr^^uiti-w turns on. Since the number density of electrons in 
matter is nominally half that of nucleons, < 1/60 ~ 2% of multi-W events 
are due to neutrino-electron interactions. 

Concerning neutrino flux attenuation due to competing processes, we 
neglect generic z/e~ weak interactions compared to generic uN weak inter- 
actions. In the energy range of interest the neutrino-electron cross section 
due to generic weak interactions is < (9(5%) of the corresponding generic 
neutrino-nucleon cross sections [BH] for z/g, z/^ and Uf^. The only exception is 
the Glashow resonance t'e + e" — > W~ at Ep^ ~ 6.3 PeV which is 0(10^ — 10^) 



times larger than the generic a^^jy [^. However, since Ue are anticipated to 



make up only (9(1/6) of ultrahigh energy neutrinos of AGN origin and the 
Glashow resonance is relevant only near the lowest of multi-W thresholds 
(before cr^^i^j.-w has fully turned on), we neglect this effect for multi-W pro- 
duction. 

We divide our discussion of neutrino-induced multi-W phenomena into 
four sections. Given that AGN are the most plausible source of our as- 
sumed neutrino flux, we investigate in sect. the conditions under which 
large neutrino cross sections for multi-W production are compatible with 
AGN neutrino production mechanisms. In sect. |4.2| we discuss constraints on 
neutrino-induced multi-W production from the Fly's Eye and in sect. [4.3| we 
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consider the possibility of using conventional air shower arrays. In sect. ^.4| we 
evaluate the potential of subsurface detectors to observe contained multi-W 
phenomena and discuss the detection of distant multi-W phenomena through 
searches for energetic muon bundles. 

4.1 AGN neutrinos with large cross sections 

Active galactic nuclei are natural candidates for ultrahigh energy neu- 
trino production In the AGN model of Stecker et al. |]T^ charged 



pions are produced in reactions such as p7 A+ nn^ when protons 
accelerated by a spherical accretion shock |^ collide with the dense gas of 



ultraviolet photons in the innermost region around the central black hole. 
Charged pions decay and give rise to neutrinos, whereas photons produced 
through vr'' decay cascade to lower energies, eventually appearing as X-rays. 
If it is assumed that the diffuse X-ray background is primarily from AGN, 
the observed X-ray flux can be used to normalize the calculation of the neu- 



trino flux. Szabo and Protheroe [|16[ have extended this model by including 
pion production through pp interactions. The AGN model by Biermann and 
collaborators fl^ differs from the model used in refs. |]15|, [l^ mainly in 



the geometry; the shocks needed for the acceleration of the protons are as- 
sumed to arise in the bipolar outflow of gas and plasma perpendicular to the 



accretion disc 43 



Among other considerations, prolific neutrino production by AGN is a 
function of the matter density in the vicinity of the pion production; if the 
medium is too dense, charged pions undergo subsequent hadronic interactions 
instead of decaying to give neutrinos. For example, for PeV pions the matter 
density should be less than 10~^ g/cm^ Assuming that the prevailing 
conditions are indeed conducive to neutrino production, we must consider 
whether or not potentially large neutrino-induced multi-W cross sections 
permit produced neutrinos to escape. For example, one might worry about 



13 



neutrino reabsorption due to multi-W phenomena through coUisions with the 
dense gas of ultraviolet photons required by the AGN model of Stecker et 
al. However, even a low parton-parton multi-W threshold of a/sq = 2.4 TeV 
corresponds to a neutrino threshold energy (colliding with ~ 40 eV photons) 
of 3 X 10^^ GeV; we can safely ignore this effect. More important is the effect 
of neutrino reabsorption by matter; the relevant parameter is the column 
density seen by particles escaping from the inner regions of AGN. 

In the Stecker et al. model the effective escape column density for neutri- 
nos from the central region of the AGN is of the order of the column density 
of the X-ray emitting region, Xfg^'j^pf ^ Xx-ray — O(10~^ — 10~^ g/cm^) 
1^, ^ whereas in the model of Biermann et al. Xf^^l^^"'^ ~ 10^ g/ cm^ 



is of the order of a hadronic interaction length [46|. As can be seen in Fig. 10 



which plots the neutrino interaction length = nip/ a'^^^^_^^{Ey), signifi- 
cant reabsorption of neutrinos by AGN is not an issue since Xy ^ Xescape- 
A related point, but one which we do not address in this paper, is the im- 
plication of large neutrino cross sections for stars near the cores of AGN. 
Neutrino interaction lengths of 0(10^ g/cm^) due to generic charged cur- 
rents for i?,^ ~ 1 — 100 PeV are sufficient to disrupt stellar evolution near 
the cores of AGN E^, ffH]; shorter neutrino interaction lengths implied by 



multi-W phenomena may have interesting consequences. 

For definiteness, we use the (revised) Stecker et al. AGN neutrino flux 



15] in the following sections for estimates of expected rates of multi-W phe- 



nomena. It should be noted, however, that the fluxes calculated in refs. 



1^, generally agree| with Stecker et al. [T^ above .1 PeV, which is the 



energy range we are interested in. In this sense our use of the Stecker et al. 
flux is intended to be representative of a large class of AGN flux models. 



2/ 



whi 
Iftsee, e 



At lower energies refs. |16, |l7|, which take pp interactions into account, give consid- 



erably larger fluxes than ref. [1^ (see, e.g., Fig. 11 in ref. [ |18| ) 
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4.2 The Fly's Eye experiment 

Independent of any neutrino flux model, the Fly's Eye array |^ puts 
upper limits on the product of the flux times total cross section for weakly 
interacting particles in the range 10^ GeV < < 10^^ GeV assuming that 



such particles initiate extensive air showers deep in the atmosphere |4^. The 
limits are deduced from the non-observation of downward-moving air showers 
within the Fly's Eye fiducial volume such that the shower axis is inclined 
80° to 90° from the zenith at the point of impact on the Earth. Showers 
meeting these criteria could only have been initiated by particles typically 
penetrating more than 3000 g/cm^ of atmosphere before interacting, which 
excludes showers initiated by ultrahigh energy photons and hadrons. 

Assuming that the weakly interacting particles referred to by the Fly's 
Eye are neutrinos, we denote the relevant cross section by a^j^ {Eu) which 
receives contributions from both multi-W and familiar charged current weak 
interactions. The Fly's Eye limits may be summarized by (ji/0"^o^)Fiy's Eye 
< 3.74 X 10-^2 X [EJI GeV)-i-^s s^^ sr-^ GeV-^ [||, ||]. Since these limits 
neglect the possibility of flux attenuation in the upper atmosphere due to 
large inelastic cross sections, they nominally apply only if cr^j^ {E^) < 10 /ib. 
If one considers a particular neutrino flux model j^"'^°^{Ey) the Fly's Eye 
limit excludes regions in the {E^^a^^) plane bounded by 



(i'^'^tcrt )Fly's Eye „vN ( \ ^ in ,,K 
< (^toAEu) < 10 /ib^ 



10^ GeV < E„ < 10" GeV. 



(8) 



If we use the AGN neutrino flux of Stecker et al. [jT5[ (i.e., set j™°*^''' = 
jStecker et a«.^^ -j-j-^g Fly's Eye cxcludcs the hatched region of Fig. ^ 

In order that {Ey) avoid the region excluded by eq. ||, only certain 



combinations of yso and o"o are consistent. For example, as shown in Fig. 11 



for ^/sQ = 8 TeV the range .5 /ib < 6"o < 81 /ib is excluded. Similar consid- 
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erations for other values of y/so result in the excluded region labeled "AGN 



u" in Figure 0. As may be quickly verified, the upper left boundary of 
the excluded region in Figure ^ corresponds to limiting situations in which 
a^^ {Ey = 10^ GeV) = 10 /ib. In principle the 10 /ib upper bound on the 
neutrino-nucleon cross section in eq. | could be enlarged by taking into ac- 
count flux attenuation in the upper atmosphere, which has been neglected 
in ref. [Q. As a consequence, one could most likely extend the excluded 
region in Fig. |12| into the upper left hand corner which, taken literally, is not 
constrained by the Fly's Eye. A further improvement of the Fly's Eye limit 
(ji'^tot^) Fly's Eye, by a factor of 10-50, is expected from the High Resolution 
(HiRes) Fly's Eye |0|. 

Active galactic nuclei are not the only conjectured sources of ultrahigh 
energy neutrinos. As shown in Fig. ^ when the proposed AGN neutrino flux 
dies off beyond ~ 1 EeV, the dominant component to the neutrino flux may 
be due to protons scattering inelastically off the 2.7 K cosmic background 
radiation (GBR) producing charged pions that subsequently decay and 
produce neutrinos |52|, |5^. The photoproduced neutrino flux, j^ '' ^, shown 
in Fig. ^ is taken from ref. |TB[. It is amusing to consider how the Fly's 
Eye constraints on neutrino-induced multi-W production are modified if we 
account for the possibility of such photoproduced neutrinos. If one takes 

^•model ^ ^-Stacker et al. ^^2.7 K g p^y^ excludes the (^^, <o^) 

region shown Fig. |13| and enlarges the excluded region in {^/so, ctq) space by 



the area labelled "2.7 K Photoproduced z/" in Fig. [12. 

Though the appearance of an enlarged excluded region is welcome, it 
is sensitive to details of the assumed GBR flux. Had we assumed a GBR 
flux component j^ '' ^ which was a factor of ten smaller than that shown in 
Fig. I (corresponding to a lower redshift), the quantity (ji/0'to^)Fiy's Eye/j™°'^^^ 
(corresponding to the solid curve in Fig. p!3| ) would not have dipped below 
10 /ib and thus would not have introduced a constraint. We should keep such 
uncertainties in mind to avoid attaching undue significance to the excluded 
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regions in Fig. |I2[ Nevertheless it is intriguing to speculate about detecting 
CBR neutrinos via multi-W processes since the prospects for detecting such 
neutrinos through generic weak interactions is poor unless the CBR neutrino 
flux is associated with a very large redshift. 

4.3 Conventional air shower arrays 



Consider a conventional air shower array (area A) which is located at an 
atmospheric depth Xq and is sensitive to showers above a threshold i^thresh- 
For showers not close to the horizon it is straightforward to show that the 
number of multi-W showers in time t occurring in the atmosphere above the 
detector is 

Number of multi-W showers = — / dE dVL -a'^^^^^_y^{E)jy{E). 

■'^thresh cost/ 

(9) 

The solid contours of Figure correspond to neutrino-induced multi-W 
events in 10'' s above a 100 km^ array which is sensitive to showers above 



100 PeV within 60° of the zenith; the AGN neutrino flux of Stecker et al. [ITS 



has been assumed. Even though the neutrino flux is smaller than the cosmic 
proton flux in the CMC model (see Fig. ^) the contours for neutrino-induced 
multi-W phenomena cover a considerably larger region in the a/Iq — ctq plane 
than the corresponding contours for proton-induced air showers (see Fig. ^). 
This is due to a more rapid growth of the multi-W cross section for uN 
scattering compared to pN scattering (Fig. |l|) and also due to a much smaller 
competing cross section from generic charged current interactions, ~ 
0{1 nb). 



The solid contours in Fig. |1^ do not account for the efficiency of an array 
to trigger on low altitude air showers. Such considerations are crucial for 
neutrino-induced phenomena since the distribution of neutrino interactions 
essentially follows the density profile of the atmosphere. In an exponential at- 
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mosphere neutrino-induced air showers may be initiated so close to the array 
that the showers do not spread out sufficiently to trigger the array. Rather 
than confine ourselves to a detailed analysis of triggering requirements, con- 
sider the following approximation. Suppose that an array does not trigger 
on showers initiated with 500 g/cw? of the detection level. This assumption 
is reasonable for vertical showers but is somewhat pessimistic for showers at 
larger zenith angles. Contours for "triggerable" neutrino-induced multi-W 
air showers follow from eq. ^ if, in the integrand, we replace Xq/ cos 6 with 
(Xo/ cos ^- 500 g/cm2). 

The dashed contours of Fig. |14| correspond to "triggerable" neutrino- 
induced multi-W showers. It is interesting to note that the contours for 
^5 — 10 events in 10'' s for a 100 km^ surface array roughly coincide with 
the lower boundary of the Fly's Eye excluded region labelled "AGN u" in 
Fig. |l^; this effect is easily understood in terms of the relevant acceptances 
and exposure times. If, in addition to the AGN neutrino flux of Stecker et 
ai, we were to assume contributions from 2.7 K photoproduced neutrinos as 
in the previous section, we would flnd contours for 0{1) event in 10"^ s for a 
100 km^ surface array which roughly coincide with the lower boundary of the 
Fly's Eye excluded region labelled "2.7 K Photoproduced z/" in Fig. |T2[ In 
other words, assuming the same neutrino flux, the sensitivity of the Fly's Eye 
to multi-W phenomena is comparable to that of a 100 km^ surface array. For 
this reason we will not discuss the characteristics of neutrino-induced multi- 



W air showers relevant to surface arrays. Mrenna |5J] has compared the 
the features of neutrino-induced air showers and generic air showers in the 
context of composite models [^] where hypothesized colored subconstituents 
of PeV neutrinos interact with typical QCD cross sections. 

Conventional air shower arrays can also search for showers close to the 
horizon. Data from the AKENO array places limits on the existence of elec- 
tromagnetic (muon-poor) horizontal air showers initiated deep in the atmo- 
sphere ||5^, 0]. Since the AKENO data applies only to shower energies lower 
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than the multi-PeV range in which we are interested |^8[, we have not inves- 
tigated its imphcations for multi-W phenomena. Exploitation of the hori- 
zontal shower limits would require consideration of distant neutrino-induced 
multi-W processes in which only prompt muons penetrate the intervening 
atmosphere and initiate electromagnetic cascades close the surface array. 

4.4 Subsurface experiments 

Detectors deep below the surface of the Earth, be they shielded by rock 
(LVD ig, MACRO [gg, Soudan-2 |31[ etc.), water (DUMAND |1T1, NES- 



TOR [0, NT-200 [0) or ice (AMANDA |T9[) offer a unique perspective 
on neutrino-induced phenomena. In this section we investigate two possi- 
ble modes for detecting neutrino-induced multi-W phenomena using subsur- 
face experiments. We first consider the prospects for observing contained 
neutrino-induced multi-W phenomena and later turn to the detection of 
muon bundles arising from neutrino interactions in the surrounding medium. 

Aside from the energy involved, contained neutrino-induced multi-W pro- 
duction would reveal its origins by its enormous multiplicity (0(400) charged 
hadrons, (9(400) photons, and a few prompt muons and electrons). Generic 
deep inelastic uN scattering and the resonant process z/g + e~ — >• W~ —>■ 
hadrons can also give contained hadron production, but only with signifi- 
cantly lower multiplicity. 

The number of neutrino-induced multi-W events occurring inside a sub- 
surface detector volume V during a time t is 

Number of multi-W events = 

t ^— <li-w / dE^ dQUE^) e-^o'i^^C^-^)/-., (lo) 
■nip J 

where X{6, 0) is the column density of material in the {6, 0) direction between 
the detector and the upper atmosphere and p is the density of the material 
in which the neutrino interaction occurs. 
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Figure |T5| shows contours for contained multi-W events in 10^ s in a 1 km^ 
volume of water at an ocean depth of 4.5 km. This arrangement approxi- 



mates the proposed SADCO acoustic array |^ which, though designed to 
use acoustic techniques to detect the resonant process z/g + e~ W~, would 
also be sensitive to multi-W phenomena which are more energetic. Acous- 
tic techniques have also been considered for AMANDA [^] and DUMAND 
[ |6l[] . The contours in Fig. ^ consider the AGN neutrino flux of Stecker et 



al. as well as the sum of the Stecker flux with the 2.7 K photoproduced 



component from Fig. ^ 

Contours for the number of contained multi-W events in DUMAND can 
be obtained from Fig. |T5| simply by scaling the appropriate volume ratio. If 
we idealize DUMAND as a 100 m x 100 m x 250 m volume under 4.5 km 
of water (neglecting the possibility that its effective acoustic volume can 
be larger than its geometrical size), the volume ratio is ~ 1/400. Due to 
a numerical coincidence, the corresponding contours for MACRO (which we 
idealize as a 77 m x 12 m x 9 m volume at a depth of 3700 hg/cm^ below the 
surface of a spherical Earth with p = 2.6 g/cm^) may be obtained by from 
Fig. by scaling by a factor of ~ 1/40000, implying MACRO'S insensitivity 
to contained events. 

Due to the enormous energies involved, one need not to concentrate on 
completely contained multi-W reactions. Of particular interest is the ability 
for subsurface detectors to detect muons which arise from energetic neutrino 
interactions up to a few kilometers away. For distant multi-W production 
the effects of producing hundreds of hadrons and photons will have died off 
well before reaching the detector but the anticipated 2-3 muons from prompt 
W decays produced with E^, ^ C(100 TeV) and ~ 0(40 GeV) propagate 
great distances. The signature of multi-W production in this case would be 
energetic muon bundles. 

The ability to detect muons from distant neutrino reactions increases a 
subsurface detector's effective neutrino target volume dramatically and is the 
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premise upon which such detectors can act as neutrino telescopes. Consid- 
erable effort has recently been directed towards the prospects of detecting 
ultrahigh energy neutrinos (most likely from AGN) using subsurface detec- 



tors |18|. Despite their limited sensitivity to such phenomena, Frejus |62 



and Soudan-2 have already placed useful observational constraints on 
AGN flux models. 

As discussed in ref. [|13], near-horizontal muon bundles in DUMAND and 
MACRO would be characteristic of neutrino-induced multi-W phenomena. 
By concentrating on large zenith angles, one can avoid the complications from 
a large background of muon bundles from generic hadronic interactions in 
the atmosphere. The number of muon bundles containing k muons detected 
during time t by a subsurface detector of length L, width W and height H is 



Number of muon bundles 
t X d cos 9 



dN, 



-H{L + W) sin^ + LW\ cos^| 
7r 



'111 



dA dt dQ 

where the quantity in square brackets is the azimuthally averaged projected 
area of the detector. The calculation of the differential flux of muon bundles 
containing k muons, dNk^/{dA dt dQ), employs the techniques of ref. [1T3 
which are summarized in the appendix. 



We present in Fig. |16| contours for muon bundles beyond zenith angles 
of 80° for MACRO and DUMAND for the AGN neutrino flux of Stecker et 



al. IjT^. Due to our assumed production of 30 W bosons, each muon bundle 



consists of approximately 3 muons. The average muon energy {E^) entering 
the detector and the average inter-muon separation (r^) depend on and 
(Tq. For example, for (a/Iq = 4 TeV, o"o = 10 nb) one expects ^ 1.5 bundles 
per 10^ s in DUMAND with {E^) ~ C(180 TeV) and (r^) = 0(2.5 m); for 
= 4 TeV, (To = 1 /ib) one expects ~ 30 bundles per 10^ s in DUMAND 
with {E^) ~ 0(70 TeV) and (r^) = 0(3.6 m). Assuming an additional 2.7 K 
photoproduced neutrino flux component at the level shown in Fig. ^ changes 
the contours of Fig. Il6| by a negligible amount. 
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It may also be possible to constrain multi-W phenomena by searching 
for non-horizontal muon bundles and thereby enlarge the accessible region in 
a/sq-o-q space. Fig. (Fig. |T^) shows contours for muon bundles for zenith 
angles between 0° and 180° for DUMAND (MACROS) for the Stecker et al. 
AGN neutrino flux. An additional 2.7 K photoproduced neutrino flux compo- 
nent at the level of fig. § changes the 1-10 event contours for DUMAND but 
has a negligible effect on the MACRO contours. Since DUMAND is specifi- 
cally not optimized for downward muons the DUMAND rates in Fig. |l^ are 
presented as a matter of completeness rather than practicality. MACRO, 
however, is sensitive to downward muons. Whereas the inter-muon sepa- 
ration expected from generic hadronic interactions high in the atmosphere 
is typically of (9(5-10 m), neutrino-induced multi-W phenomena occur pri- 
marily inside the Earth and result in much more spatially compact muon 
bundles. 

Figure |T9| compares MACRO data for pair-wise muon separation to the 
contribution expected from neutrino- induced multi-W phenomena for (a/Iq = 
2.4 TeV, (Jo = 10 /ib). The MACRO data is taken from Fig. 4 of ref. [§3 



and corresponds to muon bundles at zenith angles 9 < 60° detected by two 
supermodules operating for 2334.3 hours. The MACRO data contains con- 
tributions from muon bundles of all multiplicities; approximately half of the 
reconstructed pairs come from = 2 muon bundles. We suggest that by 
separately examining the pair-wise muon separation in bundles with fixed 
numbers of muons {e.g., n^=3) as has been done the Frejus collaboration 
|6^ , MACRO may be able to put constraints on the existence of multi-W 



phenomena. A particularly useful signature of multi-W processes in this 
respect is the energy carried by each muon. Muons arising from multi-W 



processes in Fig. IT9| would have energies of approximately 80 TeV as they 



^Preliminary rates for muon bundles in MACRO presented in ref. ||64| included only 
AGN f^— induced multi-W processes and hence are smaller than those of fig. 18 by a factor 
of 3. 
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enter the detector and may be distinguished by mechanisms such as catas- 



trophic energy loss Though some of the region in (vsq, ctq) space to 

which MACRO is sensitive is aheady excluded by the Fly's Eye (assuming 
the same AGN neutrino flux), valuable independent constraints may already 
be possible from existing MACRO data. 



5 Summary and Conclusions 

Future hadron colliders such as the proposed SSC or LHC offer the 
best prospects for observing or constraining multi-W phenomena. A naive 
measure of an experiment's sensitivity to multi-W phenomena is the size 
of the region in (a/Iq, (Tq) parameter space accessible to the experiment by 
requiring at least one multi-W event in 10^ s of operation. By this standard 
the SSC covers the most territory (see Fig. 0). For example, at the SSC one 
would expect (9(100 events / 10"^ s) if multi-W processes were characterized 
by {\/so = 20 TeV, ctq = 1 pb). Even in the relatively noisy environment of a 
high-luminosity hadron collider, the spectacular signature of > (9(30) gauge 
bosons in a single event has no conceivable Standard Model background. In 
this sense, when applied to a hadron collider, even a naive measure of the 
sensitivity to multi-W phenomena is appropriate. 

Before the commissioning of hadron supercolliders, cosmic ray physics 
suggests alternative techniques for searching for multi-W processes induced 
either by protons or neutrinos. Taken at face value, the most optimistic cos- 
mic ray constraints on multi-W phenomena come from 1) the non-observation 
of neutrino-induced air showers by Fly's Eye (see Fig. 12) which covers the 
full range of a/Iq up to 40 TeV if o"o > (9(10 — 100 nb) and 2) searches for 
horizontal muon bundles in DUMAND (see Fig. 16) whose sensitivity extends 
to (To — 0{1 nb) if ^/lo < (9(4 TeV). Additional constraints on neutrino- 
induced phenomena may be forthcoming from limits on the occurrence of 
energetic, spatially compact non-horizontal muon bundles in MACRO or, 
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more speculatively, from large energy deposits in proposed acoustic arrays. 

Understandably, none of the neutrino-based constraints presented in this 
paper are conclusive: all presume the existence of a sizeable flux of ultra- 
high energy neutrinos from AGN or neutrinos photoproduced off the 2.7 K 
cosmic background radiation. If it should happen the required neutrino flux 
is absent, no conclusions may be drawn regarding the existence of multi-W 
processes — one would then have to wait for the advent of supercolliders to 
observe or exclude multi-W phenomena. However, it is instructive to con- 
sider an intermediate scenario in which a flux of ultrahigh energy neutrinos 
is detected in the future but is found to have interactions consistent generic 
charged current processes; this too may place constraints on the existence 
of exotic phenomena. In any case, since no experiment has yet studied the 
interactions of neutrinos with energies greater than a few hundred GeV, it 
may be premature to dismiss the possibility that PeV neutrinos have novel 
interactions. 

To avoid the additional uncertainty of whether an ultrahigh energy neu- 
trino flux exists, we have also considered the possibility of multi-W phenom- 
ena induced by a cosmic protons. Unfortunately, the tradeoff for a rela- 
tively reliable proton flux is an overwhelming competing cross section due to 
generic hadronic processes. Large competing cross sections complicate mat- 
ters in two ways. First, they drastically reduce the proton flux effectively 
available for multi-W phenomena. Even by our naive measure of sensitivity 
to multi-W phenomena, a 100 km^ surface array operating for 10'' s could 
only see multi-W phenomena if 0(10 — 12 TeV) and (Tq > 0{.l — 1 /ib) 

(see, e.g., Fig. |^). Second, large competing cross sections are the source of an 
overwhelming background of generic air showers. Given that there are opti- 
mistically (9(10^ — 10^) generic showers for every shower of multi-W origin, 
we have emphasized that multi-W showers could easily be mistaken for back- 
ground fluctuations. Taking this difficulty into consideration, the true sensi- 
tivity of conventional surface arrays to proton-induced multi-W phenomena is 
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negligible. Similar conclusions apply to other techniques for proton-induced 
multi-W phenomena such as detecting downward moving underground TeV 
muons in MACRO, or searching for the anomalous longitudinal development 
of air showers with the Fly's Eye. 

The short term outlook for constraining or detecting multi-W phenomena 
in cosmic ray physics is mixed. Without making additional assumptions (such 
as assuming the existence of a large cosmic neutrino flux) one must focus on 
proton-induced processes and conclude that current and future experiments 
are effectively insensitive to multi-W phenomena over the entire range of 
parameter space where they might plausibly exist. From this viewpoint one 
must wait for terrestrial supercolliders before conclusive constraints on multi- 
W processes arc established. While this conservative scenario may very well 
be true, an exciting alternative exists. If a sizeable flux of cosmic neutrinos is 
present, not only may AMANDA, DUMAND, MACRO, NESTOR and NT- 
200 be sensitive to them through generic weak interactions, but such detectors 
may also indicate whether multi-W processes arc real or an artifact of our 
imperfect understanding of multi-TeV weak interactions. 
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7 Appendix 



In this appendix we outline our quantitative description of multi-W pro- 
cesses. In sect. |7[5 we state the assumptions and approximations used to 
model proton- induced multi-W production and in sect. ^i2| we briefly de- 
scribe our simulation of extensive air showers. Finally, in sect. [7.3| we review 
the calculation of subsurface detection rates for muon bundles from neutrino- 
induced multi-W processes. 

7.1 Proton-induced multi-W processes 

Within our working hypothesis the proton-nucleon multi-W cross section 
is given by 



pN 

multi-W 



> / dXi dX2 -—J (To 6 y/XiX2S - J Sq 



(12) 



where fi{x) is the parton distribution function corresponding to a parton 
of flavour i carrying a proton momentum fraction x. The sum extends over 
all distinct^ combinations of quarks and antiquarks (but not gluons). We 
evaluate all parton distribution functions at a scale = and employ 
the leading-order parton distributions of Tung and Morfin (fit SL) . 



In the proton-nucleon center of momentum system (c.m.s.) where the 

total energy is ^/s, we sample the parton distribution functions to generate 

the momentum fractions xi,X2 carried by the quarks participating in the 

hard interaction. In the quark-quark c.m.s. the energy of the hard subpro- 

cess is ^/I = y/xiX2S. For definiteness we assume that multi-W processes 

^Equation |l^ corrects eq. 3 of ref. Though ref. used eq. ^ for calculations, 

^nmiti-w curve in Fig. 1 of ref. is too large by approximately a factor of 2 due to 
an incorrect double-counting of contributions from unlike partons. 



26 



produce exactly 30 W bosons. A more detailed scenario should consider the 
production of Z bosons (roughly in the ratio : W~ : Z ~ 1 : 1 : 1), 
prompt Higgs bosons, prompt photons and allow for fluctuations in the total 
weak boson multiplicity. In addition, we should, strictly speaking, ensure the 
conservation of the quantum numbers carried by the quarks participating in 
the hard interaction by allowing more than W bosons in the final state of 
the hard subprocess. In the interest of simplicity we will sacrifice complete 
consistency and forgo such refinements; compared to the large number of 
gauge bosons produced we do not expect these points to play a significant 
role in our investigation of whether or not it is feasible to observe multi-W 
processes. 

We assume that in the quark-quark c.m.s. the W boson momenta are 
distributed isotropically with each W boson carrying an energy \/S/30. A 
more detailed treatment should employ 30-body relativistic phase space and 
perhaps impose dynamic assumptions such as limited pt — 0{mw) of the 
W bosons (in analogy with limited pt in QCD). However, because of the 
rapidly falling parton-parton luminosity, the quark-quark subprocess energy 
y/I tends to lie just above the multi-W threshold y/sQ. Consequently, for 
multi-W thresholds close to the kinematic limit for the production of 30 W 
bosons {i.e., a/Iq — 2.4 TeV) there is little extra energy available to have to 
worry about the precise distribution of the W bosons in momentum space. 

We employ the Monte Carlo program jetset to decay all W bosons 
and to reproduce measured W branching fractions and hadronic multiplici- 
ties. At this stage we inhibit the decays of relatively long-lived secondaries 
such as n"^, k^, p, r], k^, rj' to allow for the possibility that they may undergo 
hadronic interactions with air nuclei in the subsequent air shower. The ex- 
tensive air shower simulator described in sect. |7]^ will determine whether 
these secondaries decay or interact. 

Aside from the decays of the W bosons we must also consider those parts 
of the colliding protons which do not participate in the hard subprocess 
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— the so-called spectator fragments. Since the quarks participating in the 
subprocess carry colour, the hadronization of the spectator fragments is not, 
in general, independent of the hadronization of the subprocess system. As a 
simplification we ignore this by point and adopt the following procedure: 1) 
We treat the multi-W production subprocess as a color singlet (and hence 
hadronize it independently by decaying all 30 W bosons). 2) We replace 
the spectator fragment originating from the cosmic proton with a nucleon 
carrying the same energy. 3) We ignore the spectator fragment originating 
from the stationary "target" proton. This ansatz, especially steps 2) and 
3), is intended to embody the essential characteristics of the leading particle 
effect. As a final step we boost all W decay products and the leading nucleon 
to the Earth rest frame and inject them into the upper atmosphere to be used 
as initial conditions for an extensive air shower. 



7.2 Air shower simulation 



We use the computer program showersim ||3^ to simulate both the elec- 
tromagnetic and hadronic components of air showers generated by multi-W 
processes. The program accounts for multiple hadronic and electromagnetic 
interactions in the atmosphere and allows for the decay of unstable parti- 
cles. The interested reader is referred to the showersim documentation for 
a detailed discussion of the program's physical assumptions. We employed 
the program in its default form with few exceptions. For the electromagnetic 
components of showers we employed the elcas.5 and tail. 2 routines which 
provide a detailed evolution of photons and electrons below 200 GeV; for 
the underlying hadronic interaction model we used the 'WOO" option which 
nominally fits SPS data. 

We also used showersim to generate samples of generic air showers in- 
duced by protons and iron nuclei. We stress that though we employed show- 
ersim exclusively, we did so only for convenience since it is not clear that 
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any available shower model, especially concerning the hadronic component of 
extensive air showers, provides an accurate representation of Nature. In this 
respect, our results for multi-W air showers, generic proton or iron showers 
may not be accurate on an absolute scale but should be reliable relative to 
each other due to common simulation techniques. 



7.3 Multi-muon detection rates 

For completeness we summarize here the ingredients of our calculations 



of multi-muon detection rates. The interested reader is referred to ref. |]T3 
for additional details. 

We characterize a subsurface detector by its vertical depth D, its geomet- 
rical size (length L, width W, height H) and a muon threshold energy -Ethresh 
which is the minimum energy required of muon entering the detector in order 
that it pass completely through the detector. For our calculations we idealize 
DUMAND (with a nine string array) as a 100 m x 100 m x 250 m volume at 
an ocean depth of 4.5 km with -Ethresh = 100 GeV. For simplicity we neglect 
the effective growth of the array size with muon energy. Similarly, we idealize 
MACRO as 77 mx 12 m x 9 m volume located at depth of 3700 hg/cm^ 
below the surface of a spherical Earth of density p = 2.6 g/cm^. For MACRO 
we assume -Ethresh = 2 GeV. Since prompt muons from multi-W phenomena 
would typically arrive at DUMAND or MACRO with energies >C»(10 TeV), 
the small values of -Ethresh used above are essentially irrelevant. 

For an isotropic differential flux of cosmic neutrinos, j,^, the differential 
flux of detected events with k muons in coincidence originating from multi-W 
phenomena is given by 



dN, 



/ dEPk^{E,X)UE), (13) 

-^thresh 



dA dt dVt 

where the total column density of matter between the detector and the upper 
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atmosphere is 



X = p 



(i?e - Df cos2 d + 2DR^ - - {R^ - D) cos 9 



+ X,tm. (14) 



The first term in the column density accounts for rock/water /ice above the 
detector and the second term is the appropriate atmospheric slant depth 
using the U.S. Standard atmosphere model. 

Pkfj,{E,X) is the probability that a cosmic neutrino of energy E, initially 
separated from the detector by a column density X, gives rise to a multi-W 
event with k muons detected in coincidence, 

(15) 

where is the number of collimated muons produced in a multi-W event 
and Na = rn'^. We assume = nw/9 — 3, corresponding to prompt muons 
from the decay of 30 W bosons. 

Pfj^{E, Ethresh, X') is the detection probability for a typical prompt muon 
produced in a multi-W process induced by a neutrino with energy E. Under 
the assumption that the ?t,vk weak gauge bosons are distributed isotropically 
in the subprocess rest frame, 

~ / iTi ZJi v'\ 1 -^/j (-^thresh 5 ) 

p^{E, Ethresh,X ) = 1 - nw— , (16) 

where E'^{Ethresh, X') is the solution of the muon energy-range relation 

X' = (17) 

yi?„«(^)+/3(^)^' ^ ^ 

We neglect the stochastic effects of range straggling which may become im- 
portant for muon energies above 10^ GeV and use parametrizations of a{E) 
and P{E) from refs. |70|, respectively. 

The distribution of the inter-muon separation in muon bundles from 
neutrino-induced multi-W processes is obtained as a byproduct of using stan- 
dard Monte Carlo techniques to evaluate the integral of eq. O which gives 
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the number of underground muon bundles. For each point contributing to 



the integral of eq. |Tl| one generates a muon bundle configuration, propa- 
gates it to the detector and calculates its contribution to the distribution 



of pair- wise separations using the integrand of eq. |Tl| as a weight. Muon 
bundle configurations are generated according to the assumptions outlined 
in sect. [7.1| . Namely, for a cosmic neutrino of energy initiating a multi-W 
process, one samples the parton distribution functions of the target nucleon 
to determine the quark-neutrino c.m.s. energy ^/I = J2mpEj^x (> -\/so)- 



Since prompt muons from W boson decay are distributed isotropically in the 
quark- neutrino c.m.s., one generates momentum vectors for three prompt 
muons in that frame and then boosts the results to the Earth rest frame. Af- 
ter including the effects of muon energy loss and multiple Coulomb scattering 
in the medium surrounding the detector, the 2 < < 3 muons reaching the 
detector determine n^(n^ — 1)/2 pairwise separations which contribute to the 



muon separation distribution such as shown in fig. [19 . 
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Figure Captions 



1. Universal curves parametrizing multi-W production cross sections in 
proton-nucleon (pN), proton-electron {pe~), neutrino-nucleon {i^N) 
and neutrino-electron (i^e~) collisions. Curves are for protons and neu- 
trinos with laboratory energy E colliding with nucleons and electrons 
at rest. 

^{pN thresh) ^ so/{2mp) IS the proton threshold energy for pN 
multi-W processes. The (pN) curve corrects an error in ref. the 
corresponding curve of ref. is too large by approximately a factor 
of 2 (see footnote in appendix). 

2. Contours corresponding to 1 and 10 multi-W events in one year (10^ s) 
of operation for the LHC {C = lO^^ cm'^ s"!) and the SSC (£ = 

^q33 ^^-2 

3. Event number contours in 10^ s for proton-induced multi-W air showers 
assuming the Constant Mass Composition model for proton flux. Solid: 
100 km^ conventional surface array sensitive to -Eshowcr > 1 PeV at 
zenith angles 9 < 60°. Dashed: Fly's Eye array sensitive to -Eghower > 
100 PeV using aperture of ref. . 



4. Differential flux of protons and neutrinos used in text. The Constant 



Mass Composition proton flux is from ref. . The diffuse neutrino flux 
from active galactic nuclei (AGN) and the 2.7 K photoproduced neu- 
trino flux are taken from ref. [|T5|. Neutrino fluxes shown are summed 



over species in the proportion z/^ : z/^ : i/g : J^e = 2 : 2 : 1 : 1. 

5. Differential flux of proton-induced multi-W air showers assuming the 
Constant Mass Composition (CMC) model proton flux for flxed multi- 
W production parameters. 

6. Lateral distributions of electrons {Ee > 1 MeV), muons (i?^ > 1 GeV) 
and hadrons (-Ehad > 1 GeV) in 5 PeV vertical air showers (lower 
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three curves in each plot) and 30 PeV vertical air showers (upper three 
curves in each plot) at atmospheric depth of 800 g/cm^. Solid curves 
correspond to proton-induced multi-W showers assuming a/Iq = 2.4 
TeV and any value of ctq. Dashed (dot-dashed) curves correspond to 
generic showers initiated by proton (iron) primaries. Each curve is an 
average over 25-100 showers including variations in the depth of first 
interaction. 

7. Average fraction of total shower energy Ep carried by generic compo- 
nent of a proton- induced multi-W air shower. E^"^"^^^ = So/{2mp) is the 
proton threshold energy for multi-W production. 

8. Lateral distributions of muons with > 1.5 TeV at atmospheric depth 
of 800 g/cm^ for vertical air showers initiated by 30 PeV primaries. The 
multi-W processes assume a/Iq = 2.4 TeV. Each curve is an average 
over 100-500 showers including variations in the depth of first interac- 
tion. 

9. Longitudinal development curves for 150 PeV hadron- induced vertical 
air showers. The depths of first interaction are fixed at 42 g/cm^ for p- 
induced multi-W showers, 42 g/cm^ for generic p-initiated showers and 
11 g/cm^ for generic Fe-initiated showers. A multi-W parton-parton 
threshold of a/s^ = 5 TeV is assumed. A number of curves are shown 
to illustrate the similarity between multi-W showers and fluctuations 
in generic showers. 

10. Neutrino interaction length due to combined effects of generic charged 
current interactions and multi-W processes. 

11. Excluded region (hatched) in E^, — space from combination (solid 
line) of Fly's Eye limits with the AGN neutrino flux of Stecker et al. [p^ j. 
Dashed lines indicate limiting cases of for (a/Iq = 8 TeV, oq = 
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.5 /ib) and (v^ = 8 TeV, ctq = 81 /wb) which are consistent with the 
Fly's Eye hmits. 

12. Regions of muhi-W parameter space excluded by the Fly's Eye. The 
region labelled "AGN z/" is excluded if one assumes only the AGN neu- 



trino flux of Stecker et al. [^. The region labelled "2.7 Photoproduced 
z/" is excluded in addition if one includes the neutrino flux contributions 
due to the cosmic background radiation shown in Fig. ^. 

13. Excluded regions in — a'^^ space from combination (solid line) of 
Fly's Eye limits with the AGN neutrino flux of Stecker et al. |]l5l 
(hatched) and the flux due to the 2.7 K cosmic background radiation 
(double hatched) (c.f. Fig. H). Dashed lines indicate limiting cases of 

for (v/l^ = 8 TeV, (Tq = 48 nb) and {^/Fo = 8 TeV, ao = 81 /xb) 
which are consistent with the limits. 

14. Event number contours for neutrino-induced multi-W extensive air 
showers (-Eshower > 100 PeV, zenith angle < 60° ) in 10^ s for a 100 km^ 
conventional surface array (vertical depth 1000 g/cm^). Solid contours 
includes all showers. Dashed contours include only showers initiated a 
minimum of 500 g/cm^ away from array. The AGN neutrino flux of 



Stecker et al. 15 is assumed. 



15. Contours for neutrino- induced contained events in 1 km'^ volume of 
water at an ocean depth of 4.5 km in 10^ s (approximately the arrange- 
ment of the proposed SADCO acoustic detector). The neutrino flux of 
Stecker et al. [l^ is assumed (see Fig. ^). 



16. Contours for neutrino-induced multi-W muon bundles at zenith angles 
9 > 80° in 10^ s at MACRO and DUMAND assuming the AGN neutrino 



flux of Stecker et al. 15 
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17. Contours for neutrino-induced multi-W muon bundles for all zenith 
angles in 10'' s at DUMAND assuming the AGN neutrino flux of Stecker 
et al. PI (see Fig. |). 



18. Contours for neutrino-induced multi-W muon bundles for all zenith 
angles in 10^ s at MACRO assuming the AGN neutrino flux of Stecker 



et al. 15 



19. MACRO pairwise muon separation data compared with expecta- 
tions for neutrino-induced multi-W phenomena for (^/sq = 2.4 TeV , 
(To = 10 fih) assuming the AGN neutrino flux of Stecker et al. [|l^ . 
MACRO data corresponds to two supermodules operating for 2334.3 
hours sensitive to bundles with zenith angle 6 < 60°. 
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Table 1 Average particle multiplicity at atmospheric depth of 800 g/cm^ 
for vertical extensive air showers generated by 5 PeV (and 30 PeV) 
primaries. Multi-W showers are proton-induced. Showers labelled p 
and Fe contain only generic interactions. 



Shower 
Type 


> 1 MeV 


> 1 GeV 


(iVhad) 

^had > 1 GeV 


E^ > 1.5 TeV 


multi-W 


1.1 X 10^ 
(1.4 X 10^) 


5.9 X 10^ 
(1.9 X 10^) 


2.2 X lO'^ 
(1.4 X 10^) 


18 

(40) 


P 


2.4 X 10^ 
(1.8 X 10^) 


4.5 X 10^ 
(1.5 X 10^) 


2.9 X lO-'^ 
(1.5 X 10^) 


5 
(20) 


Fe 


1.3 X 10^ 
(1.1 X 10^) 


8.7 X 10^ 
(3.3 X 10^) 


3.3 X 10^ 
(1.4 X 10^) 


14 
(60) 
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